Braid Group and Topological Phase Transition in Nonequilibrium Stochastic Dynamics 
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We show that distinct topological phases of a band structure of a non-Hermitian Hamiltonian can be classified 
with elements of the braid group. As an example we consider a non-Hermitian evolution of the statistics of 
stochastic transports. We show that topologically nontrivial phases have detectable properties, including the 
emergence of decaying oscillations of certain statistical characteristics, and discontinuities in the steady state 
statistics of currents. 
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Numerous physical systems are described by Hermitian 
Hamiltonians, H(x), which periodically depend on a vector 
of parameters x with a period 2tt. As functions of %, the 
eigenvalues A^ (%) of H(x) represent the band structure of the 
system. Well known examples include Bloch bands for elec- 
trons and phonons in periodic crystal potentials [1] or Floquet 
bands for periodic-driven quantum dynamical systems [2, 3]. 
If we trace the dependence of initially nondegenerate eigen- 
values of a Hermitian Hamiltonian on %, it is possible that, at 
some %, a pair of bands Xi(x) an d Xj (%) goes through an ex- 
act crossing point. Such events correspond to quantum phase 
transitions [4] that indicate, for example, a change of topolog- 
ical order of Floquet bands in periodic momentum space [3], 
or electronic Bloch bands in periodic position space [5, 6] and 
their phononic/magnonic counterparts [7]. 

Dissipative open systems described by non-Hermitian 
H(x) a l so possess the band structures [8, 9], once they have 
periodic dependence on position, momentum, or others. In- 
vestigation of these topological band structures is currently 
an active research field [10-12]. In this Letter, we demon- 
strate that topological properties of the band structure of a 
non-Hermitian Hamiltonian can be classified with elements 
of the braid group (braid diagrams). In topologically nontriv- 
ial phases such systems can undergo a spontaneous symmetry 
breaking and show specific detectable effects. 

Consider first a Hermitian Hamiltonian with a periodic- 
ity H(x) = H(x + 27r) along one of the parameters 
and with nondegenerate eigen- spectrum Xi(x), i = 1 • . .N. 
Such real eigenvalues of H(x) must have the 27r-periodicity: 
Xi(x) = Xi(x + 27r). In this Letter, we will explore con- 
sequences of the observation that in the non-Hermitian case 
the latter restriction generally does not hold: Eigenvalues of 
a non-Hermitian Hamiltonian are generally complex. Hence, 
the periodicity of H(x) oni y implies the periodicity of the 
whole nondegenerate complex eigen- spectrum but does not 
mean that a single eigenvalue has such a periodicity. For ex- 
ample, one may encounter a situation when 

A,(x + 27r) = A i (x), A i (x+27r) = A,( X ), (i^j), (D 

i.e. two complex eigenvalues, A$ and Aj, twist in a braid so 
that the 2tt periodicity for an individual eigenvalue is lost. 
This situation cannot be met by nondegenerate real Xi(x) be- 
cause such a case would inevitably lead to exact crossing of 
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FIG. 1. (color online). Two band structures. In order to satisfy Eqs. 
(1), bands should either (a) cross or (b) create a braid pattern. The 
latter is possible only when eigenvalues are allowed to be complex 
(for non-Hermitian case). 



two eigenvalues at some as illustrated in Fig. 1(a). In con- 
trast, Fig. 1(b) demonstrates that complex eigenvalues can 
avoid each other when x changes continuously while satis- 
fying Eqs. (1). It is straightforward to generalize this ob- 
servation, namely, the twisting pattern of N band of a non- 
Hermitian N x N Hamiltonian as functions of the parameter 
X can be generally mapped to some braid diagram of the braid 
group B N . 

Two band structures are topologically equivalent if all their 
N bands can be continuously deformed into each other with- 
out encountering an exact eigenvalue crossing. Different el- 
ements of the braid group cannot be continuously deformed 
into each other. This means that topologically equivalent 
bands are described by the same diagram of the braid group, 
i.e. the braid group provides a natural classification of dif- 
ferent patterns of the band structure in systems with non- 
Hermitian Hamiltonians. Varying parameters of the Hamil- 
tonian, the topology of the band structure can only change 
at a crossing of a degeneracy of eigenvalues. We will call 
such topological phase transitions the braid phase transition 
(BPT). 

The fact that a specific band in a topologically nontriv- 
ial phase may not possess the periodicity of the Hamiltonian 
provides the means for the spontaneous symmetry breaking 
(SSB) phenomenon [13], i.e., given a symmetry of the equa- 
tions of motion, a physical system can prefer certain solutions 
that are not invariant under this symmetry. The full symmetry 
is preserved only in the ensemble of all solutions. The SSB 
is a general principle that underlies a vast number of physical 
phenomena, ranging from ferromagnets and superconductors 
to the origin of masses in elementary particle physics [14]. We 
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will show that BPTs can lead to an unusual mechanism of the 
phase transition that is in some sense opposite to SSB, namely, 
while individual bands in a nontrivial phase do not respect the 
symmetry of the Hamiltonian, the observable characteristics, 
such as the moment generating function of currents, may keep 
the original symmetry. In such a case, the system may, for 
physical reasons, choose different solutions at different values 
of the periodic parameter x » which results in discontinuity of 
such observables as functions of x- We will also show that, 
similar to Goldstone modes in SSB phases, nontrivial braid 
phases correspond to the appearance of certain modes with 
new properties. 

As the proof of principle, we will explore the band struc- 
ture of the operator H(x) that describes the transport statis- 
tics in two elementary kinetic models. The so-called counting 
parameter x conjugates to the discrete transport quantities of 
interest, e.g., the particle number, and thus leads to the 2tt 
periodicity of H(x)> which is just similar to the periodic Bril- 
louin zone resulting from the discrete lattice spacing in solid 
state physics. 

We first consider a Markovian stochastic transport of par- 
ticles between two reservoirs through an intermediate system 
that we will call a bin [see Fig. 2(a)]. When the bin is empty, 
a particle can jump from the left (right) reservoir into the bin 
with kinetic rate fa{k±). When the bin has a particle inside 
it, the particle in the bin escapes to the left (right) reservoir 
with rate ks (fa)- This model is known in biochemistry as 
the Michaeils-Menten model [15, 16]. Statistics of events in a 
system that corresponds to this model was measured in [17]. 
Similar stochastic models also appear in description of statis- 
tics of transports through nanoscale electric circuits [18] or 
molecular phononic devices [19]. 
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FIG. 2. (color online), (a) The Michaelis-Menten-like model of a 
bin coupled to two particle reservoirs. The wavy line indicates that 
transitions related to the right reservoir are monitored, (b) Braid 
phases of the counting statistics behavior of the kinetic model at 
k± = &2 = 1. The dashed line indicates the detailed balance condi- 
tion fcife = fe&4. 

Let us introduce two sets of joint probabilities, po(n,t) 
and pi(n,t) of that during time interval [0,t], there will be 
exactly n particles transferred into the right reservoir while 
at time t the bin is, respectively, empty or filled with a 
particle. The master equations for Po/i(n,t) are given by 
Po(n,t) = -(fa +k4)po(n,t)-\-k 2 pi(n-l,t) + kspi(n,i), 
Pi(n,t) = -(k2 + ka)pi(n,t) + kipo(n,t) + k4Po(n + l,t)- 
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FIG. 3. (color online). Braid diagrams of band structures for different 
phase regions in Fig. 2. (a) In phase I (fe = 1.5, k± = 0), both 
eigenvalues are periodic functions of % with period 2tt. (b) In phase 

II (&3 = 0.5, &4 = 0), eigenvalues never cross but twist to satisfy the 
condition Ai(x + 27r) = A2(x) andA2(x + 27r) = Ai(x). (c) Phase 

III (fe = 1.5, /c4 = 1.5) is similar to phase II, but with a reversed 
chirality. 

Upon a Fourier Transform, they translate into -^(zq^zi) t = 
#(x)(>o,£i) T where z /i(x,t) = Er=-oo^o/i(M)e ixn 
with the effective Hamiltonian H(x) given by 

tt( \ ( -fa - fa fae %x + fa\ 

To describe the full counting statistics of currents into the 
right reservoir, the moment generating function, Z(x,t) = 
z o(Xj t) + z i(X: t)> can be obtained as 

Z(x,t) = (l|exp(ff( X )t)|p>, 0) 

where (1| = (1, 1), and \p) = (po,Pi) T is the vector of the 
initial probabilities. The eigenvalues of H(x) are given by 

A±(x) = g ' (4) 

where K = k ± + fa + fa + fa, e± = e ±ix - 1. BPTs 
correspond to exact crossing of two eigenvalues, which can be 
found by equating the expression in the square root in (4) to 
zero. This corresponds to the condition, x = tt, and constraint 
on kinetic rates: (k\ + fa + fa + k^) 2 — Sfafa — 8fak^ = 0. 
This equality describes the phase boundaries, which divide 
the space of positive ks and k± into several phase regions, as 
shown in Fig. 2(b). 

In order to show that the phase transition corresponds to the 
topology change of the braid diagram of the band structure, 
we plot eigenvalues (4) as a function of x £ [0, 2tt) in Fig. 3. 
Fig. 3(a) corresponds to a choice of parameters in phase I, with 
a trivial periodicity. Fig. 3(b) corresponds to phase II with 
eigenvalues that are not periodic but rather twisting around 
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FIG. 4. (color online). Manifestation of BPT by the different time 
evolving patterns of the parity probability Z(ir) in different phases. 
Phase I (&3 = 1.5, = 0) shows monotonic decay, while phase 

II (fe = A)4 = 0) exhibits oscillatory decay. To expose oscillations 
of phase II, Z(ir) was multiplied by exp(Kt/2) in the inset. Phase 

III shows the similar behavior as phase II. The initial condition is 

|p> = (o,i) T 



into each other during a continuous change of \ from to 2tt 
and satisfying Eqs. (1). Fig. 3(c) corresponds to phase III, 
whose braid has a reversed chirality compared to phase II. 

The steady state current of the stochastic transport is given 
by di x \+(x)\ x =o, from which we see that different chirali- 
ties of the braid diagram correspond to different directions of 
transports. In phase II in Fig. 3(b), the particle current flows 
from left to right while in phase III in Fig. 3(c), the particle 
current flows from right to left as a consequence of the re- 
versed chirality of the braid diagram. The detailed balance 
case, fci&2 = ksk^, always belongs to the topologically trivial 
phase. It describes the equilibrium condition with zero cur- 
rent, which splits the phase diagram in Fig. 2(b) into two parts 
with opposite current directions. 

Across the phase boundaries between phase I, II, and III, the 
topology of the braid diagram will experience abrupt change 
via transitions through the degeneracy of bands. The fact that 
the degeneracy point is encountered at a finite value of the 
counting parameter, x = ^, means that it influences properties 
of the full counting statistics rather than the lowest cumulants 
of a current distribution. In fact, at \ — n me matrix H(tt) 
has all real entries and the corresponding moment generating 
function (We will call it the parity probability) 

CO 

Z(Tr,t)= [Pjpn,t)- Pj (2n+l,tj\, (5) 

n— — oo j =0,1 

has a simple physical meaning [20], namely, it is the differ- 
ence of the probabilities to observe an even and an odd num- 
ber of transitions into the right reservoir during [0, t\. 

The BPT directly influences the time-dependent behavior of 
Z(tt): Let (v± | and \u±) be, respectively, the left and the right 
eigenstates of the matrix H(tt), corresponding to eigenvalues 
e± = A±(tt). Eq. (3) gives then Z(tt) = C+e £ + t + C_e £ - t , 
where C± = (l\u±)(v±\p). Since both eigenvalues e± have 
negative real parts, the absolute value of Z(ir) decays expo- 
nentially with time. In phase I, both eigenvalues are pure real 
and s+ > £-. Hence the term C-e £ ~ l quickly becomes ex- 
ponentially suppressed in comparison with C+e £+t . At time 




„ x 0.0 0.5 1.0 1.5 2.0 2.5 3.0 

(a) (b) k x 

FIG. 5. (color online), (a) Scheme of three-state cyclic kinetics. The 
wavy line indicates that transitions between state 1 and 2 are moni- 
tored, (b) The phase diagram of the transport statistics of the three- 
state model. The dashed line indicates the detailed balance condition 
&ite&3 = k^k^ks at equilibrium. 



scales larger than 1 / (e + — S- ) the behavior of Z(ir) is then to- 
tally described by a single exponent, as we illustrate in Fig. 4. 

Consider next the case of parameters in phase II. Two eigen- 
values, e±, are now complex conjugate to each other. They 
have the same real part, which means that it is impossible to 
disregard one solution in favor of another one. Moreover, it is 
easy to see that C_ = C+ and the time-evolution of the parity 
probability becomes oscillating 



Z(tt) = Ae- Kt/2 cos(Im(e+)t + O ), 



(6) 



with some constants A and <po that depend on initial condi- 
tions. In Fig. 4, we show numerical results of the evolution 
of the parity probability in phase II, which confirm that the 
BPT to the topologically nontrivial phase corresponds to os- 
cillations of Z(tt). This also implies that the conventional 
assumption of sufficiency of only one eigenvalue to estimate 
the full counting statistics at long times cannot be used in the 
nontrivial braid phase II. Phase III shares similar behaviors. 

As a more complex example, we consider next a Markovian 
three-state kinetic model shown in Fig. 5(a). It is widely used 
as a model for cyclic enzyme reactions [21] and molecular 
motors [22]. With the transitions between state 1 and 2 being 
monitored, the evolution of the counting statistics is described 
by ^(z u z 2j z 3 ) T = H(x)(z u z 2j z 3 ) T , where the effective 
Hamiltonian H(\) is 



H(X) 




(7) 



and = EZ-ovPji^ty* 71 ' U = 1,2,3); Pj (n,t) 

is the joint probability for a system to be in state j while there 
are already n transitions from state 1 to 2 observed by time t. 

The system now has three bands Aj=i ? 2,3(x) an d a richer 
braid phase diagram. For the sake of simplicity, we fix 
ks = 2.1, &4 = 2, &5 = 0.1, &6 = 0-1 an d only plot the phase 
diagram in the space of (fci, k 2 ) in Fig. 5(b). There are four 
braid phase regions [23]. Phase I is the trivial braid with each 
band keeping the 2tt period [see Fig. 6(a)]. Phase II can be de- 
veloped from phase I by crossing Ai (x) and A2(x) at X = ^. 
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FIG. 6. (color online). Braid diagrams of band structures for different 
phases of three-state cyclic kinetics in Fig. 5. (a) Phase I with fa = 
2, k 2 = 0.1. (b) Phase II with fa = 1, k 2 = 3. (c) Phase III with 
fa = \ l k 2 = 1. (d) Phase IV with fa = 2, k 2 = 2. 



The two new bands twist and transfer into each other after a 
continuous change of x by 2n, with the third band keeping its 
periodicity intact [see Fig. 6(b)]. The nontrivial topology of 
bands is manifested by the decaying oscillation of the parity 
probability Z(ir) = Y?j=i z j ( k k )- 

In phase III, the situation is different. The first band re- 
mains intact while the second and third bands are twisting so 
that A 2 (x + 2tt) = A 3 (x) and A 3 (x + 2tt) = A 2 (x) [see Fig. 
6(c)] . This braid phase can be reached from phase I via the de- 
generacy of A 2 (x) an d A 3 (x) at x = 0, which then splits into 
a conjugate eigenvalue pair. Since at x = 0, Zj (x = 0) de- 
notes the conventional probability of the system to be in state 
j, this new topological phase corresponds to the emergence 
of decaying oscillations of conventional state probabilities to 
their steady states. 

Phase IV is a combination of phase II and III, in which three 
bands are twisted around each other. Each band becomes pe- 
riodic in x with a period equal to 6tt. The observable effect 
of this braid phase is also a combination of decaying oscilla- 
tions of both the conventional state probabilities and the parity 
probability Z(ir). 

Now, we would like to discuss the non-analytic behavior of 
the moment generating function near x = 7r, which has also 
encountered in other contexts previously [20, 24-27]. In par- 
ticular, in quantum ID electron systems it was observed that 
the analytical continuations of generating functions at x = 
and 2tt do not coincide near % = 7r [25]. It was proved re- 
cently [27] that near x = tt the generating function should be 
given by a sum of the analytical continuations of both solu- 
tions. In our model, this fact has a simple explanation. The 
parity probability in phase II is, indeed, given by a sum of two 
equal in absolute value contributions because two eigenvalues 
of H(tt) have equal real part in this phase. In phase II, an in- 
dividual eigenvalue near x = cannot be uniquely continued 
to the point X — ^ because the square root in (4) is doubly 
valued. At the steady state and at x 7^ tt, the system chooses 
the solution with the largest real part of the eigenvalue, which 
is the part of one band at x £ [0, 7r) and the part of another 
band at x £ (tt, 2tt], Hence, the appearance of the disconti- 



nuity at x = is the result of the selection of the eigenvalue 
with the largest real part at long observation time limit. 

In equilibrium statistical physics, a phase transition can 
be identified by the non-analytic singularity of free energy, 
e.g., Lee- Yang zeros [28]. In the nonequilibrium transports, 
the moment generating function Z(x,t) of dynamical observ- 
ables is the counterpart of the partition function and lnZ(x) 
is the counterpart of the free energy. Thus, the zeros in Z(x) 
and their time-dependence play a similar role to the Lee- Yang 
zeros and can be used to explore the dynamical phase transi- 
tions [29-33]. Indeed, as shown by Eq. (6), in the nontrivial 
braid phase region the parity probability Z(tt) oscillates and 
its zeros as functions of time are one-to-one mapped to the 
Lee- Yang zeros of the stochastic transport system. 

Finally, it is worth remarking that out- of- equilibrium is 
the necessary condition for the BPT because when de- 
tailed balance is respected, the eigen-spectra of the non- 
Hermitian Hamiltonian with counting parameters are always 
real. At equilibrium, due to the detailed balance condi- 
tion, the nonequilibrium force (affinity) A is zero (in our 
two examples, A — ln[jj^|j] for Michaelis-Menten-like ki- 
netics and A = ln[ ^^^ ] for three-state cyclic kinetics). 
Therefore, the Gallavotti-Cohen symmetry of the spectrum 
^j(x) = ^j(~X + 2*4.) [34] reduces to the even symmetry 
^j(x) = x)- Consequently, the odd terms of %x in the 
eigen-spectra disappear and only even terms survive so that 
the spectra are all real [35, 36] and band structures can not 
form nontrivial braids. To have nontrivial braids, the detailed 
balance condition at equilibrium must be broken. 

In summary, we have shown that, when parameters are 
varying, the band structures of non-Hermitian Hamiltonians 
may encounter topological phase transitions that can be clas- 
sified by the elements of the braid group. In application to the 
full counting statistics of currents, we found that topologically 
nontrivial phases correspond to oscillating behavior of certain 
statistical characteristics that is detectable by today's measur- 
ing technique [37]. They also correspond to discontinuities 
of the moment generating function of currents induced via se- 
lection of different band branches in the limit of long time 
evolution. Open problems include generalizations of BPTs to 
Hamiltonians that depend on a multi-component Bloch vector 
and the manifestations of BPTs in quantum mechanical full 
counting statistics and Floquet theory. 
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